Thermal entanglement of Bosonic atoms in an optical lattices with nonlinear couplings 
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The thermal entanglement of two spin-1 atoms with nonlinear couplings in an optical lattices is 
investigated in this paper. It is found that the nonlinear couplings favor the thermal entanglement 
creating. The dependence of the thermal entanglement in this system on the linear coupling, the 
nonlinear coupling, the magnetic field and temperature is also presented. The results show that the 
nonlinear couplings really change the feature of the thermal entanglement in the system, increasing 
the nonlinear coupling constant increases the critical magnetic field and the threshold temperature. 
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Entanglement as a valuable resource for quantum in- 
formation processing (QIP) 0, 0, has attracted a lot 
of attention in recent years, from both experimental and 
theoretical studies 0|. Since the entanglement is fragile, 
the problem of how to create stable entanglement remains 
a main focus of recent studies in the field of quantum in- 
formation processing. The thermal entanglement, which 
differs from the other kind of entanglement by its advan- 
tages of stability, requires neither measurement nor con- 
trolled switching of interactions in the preparation pro- 
cess, hence the thermal entanglement in various systems 
is an attractive topic and worth intensively studying. 

The system of atoms in optical lattices is among the 
promising candidates for quantum information process- 
ing. It may take the advantage of the technology used in 
atom optics and laser cooling based on the optical ma- 
nipulation of atoms []|. Besides, it also holds the merit of 
eventual possibility to scale, parallelize and miniaturize 
the device in QIP. 

The thermal entanglement has been extensively stud- 
ied for various systems including isotropic Heisenberg 
chain d, 0, E|, anisotropic Hensenberg chain Ising 
model in an arbitrarily directed magnetic field |10| . and 
cavity-QED [ill since the seminal works by Arnesen et al. 
and Nielsen |l2j| . Based on the tools developed within 
the context of quantum information theory, the relax- 
ation of a quantum s yste m towards the thermal equilib- 
rium is investigated [13j, which provides us a different 
mechanism to model a system arriving at the thermal 
entangled states. For a specific Heisenberg chain in con- 
densed matter physics, the only ranging variables are the 
magnetic field and the temperature jfj, in cavity-QED 
system, the exchange constant (the linear coupling in 
our case) is adjustable in addition to the temperature 
and magnetic field. The development of laser cooling 
and trapping provides us more ways to control the atoms 
in traps. Indeed, we can manipulate the atom-atom cou- 
pling constants and the atom number in each lattice well 
with a very well accuracy [lj, [l5| . 

In this paper, we study the thermal entanglement in 
optical lattice with nonlinear couplings. We calculate the 



thermal entanglement as a function of the nonlinear cou- 
pling constant, linear coupling constant, the temperature 
as well as the external magnetic field. We will confine 
ourself in this paper to the case of K < and J < that 
is relevant to the recent experiment conducted on 23 Na 
atoms. As we will show you later on, there is no ther- 
mal entanglement in the regime of K > J when J < 
similar to the results for the isotropic Heisenberg model 
[a] . Our studies also show that the critical magnetic field 
and the threshold temperature is obviously increased by 
the presence of the nonlinear couplings. 

Our system consists of two wells in the optical lattice 
with one spin-1 atom in each well. The lattice may be 
formed by three orthogonal laser beam, and we may use 
an effective Hamiltonian of the Bose-Hubbard form 0] 
to describe the system. The atoms in the Mott regime 
make sure that each well contains only one atom. For fi- 
nite but small hopping term t, we can expand the Hamil- 
tonian into powers of t and get (lij . 

H = e + J(S 1 -S 2 )+K(S 1 -S 2 ) 2 , (1) 

where J = —^-,K = — — jf- with t the hopping 
matrix elements, and e = J — K . U s (s = 0, 2) represents 
the Hubbard repulsion potential with total spin s, a po- 
tential U s with s = 1 is not allowed due to the identity 
of the bosons with one orbital state per well, Si(i = 1,2) 
denote the spin vector Si = (Six, Si Vj Si Z ). This Hamil- 
tonian differs from the usual Heisenberg model by the 
nonlinear couplings. Since term e contains no interac- 
tion, we can ignore it in the following discussions and 
it would not change the thermal entanglement. In the 
presence of external magnetic field, the Hamiltonian Eq. 
Q becomes 

H = J(Sx ■ S 2 ) + K{S! ■ S 2 ) 2 + B(S U + S 2Z ), (2) 

where the magnetic field B along the z-direction is as- 
sumed. When the total spin for each site Si = l(i = 1, 2), 
its components take the form 
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The ground state of Hamiltonian Eq.(J5J is expected to 
be the dimer phases, as the most recently study 01 con- 
clude, this is quite different from the Heisenberg chain 
without nonlinear couplings. And the nonlinear cou- 
plings would make the thermal entanglement different 
from that in the usual Heisenberg model, too, as you will 
see. 

To get the thermal entanglement, we first present 
the eigenvalues and the corresponding eigenstates of the 
Hamiltonian Eq.JSJ), 
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(|0,-1) + 1-1,0)), 
(|1,0) + |0,1)), 
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E 6 = Jf-J;|* 6 ) = -L(|l,-l)-|-l,l)), 

E 7 = K + J + 2B;\^ 7 ) = |1,1), 

E 8 = K + J-2B;\V R ) = |-1,-1), 

E 9 = 4if-2J;|* 9 ) = -±=(|l,-l) + |-l,l)-|0,0)). 

(4) 



The state of the above system at thermal equilibrium 
is p = Z^ 1 exp(— f3H), where Z = Tr[exp(— f3H] is the 
partition function and (3 = 1 /fc^T (ks is Boltzmann's 
constant. For simplicity we will set fcg = 1 hereafter). 
In terms of the eigenstates and the corresponding eigen- 
values, the state of the system can be expressed as 
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with the partition function 



Z = 2e- pK cosh(f3B)(l + 2 cosh(/3J)) + 2e' f)(K+J) cosh(2/?5) + e ~^ K ~^) . 
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We will choose the negativity as the entanglement mea- 
sure [HI, 

V Ta II - 1 
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where lip 7 '' 4 || denotes the trace norm of the partial trans- 
pose p A . The negativity N(p) is equivalent to the ab- 
solute value of the sum of the negative eigenvalues of 
P Ta 0] . Although the negativity lacks a direct physical 
interpretation, it bounds two relevant quantities in quan- 
tum information processing-the channel capacity and the 
distillable entanglement. As the negativity is a com- 
putable measure of entanglement for bipartite system 
with any dimension, we here choose it to measure the 
thermal entanglement. 

We have performed extensive numerical calculations 
for the entanglement measure, some selected results are 
presented in figures from 1 to 3. Figure 1 shows the plot 
of the negativity as a function of the magnetic field B and 
temperature T. For B — 0, the state j^g) is the ground 



state an d the others in Eq. (0J are excite states. In 
this case, the maximal entanglement is at T = and it 
decreases with T due to mixing of the excited states with 
the ground state. For a higher value of B greater than 
B c = 3/2(J—K), I^Pg) becomes the ground state. In that 
case, there is no entanglement at T = 0, but we may in- 
crease the entanglement by increasing T, that is to bring 
entangled eigenstates such as \*f?i) (i = 1,2,3,4,5,6,9) 
into mixing with the ground state. It is interesting to 
note that the critical field B c depends on both the linear 
coupling J and the nonlinear coupling K. With K = 0, 
B c gives rise to 3/2 J that means no entanglement in this 
case (i.e. for J < 0) at any temperature and with any 
values of magnetic field. We would like to address that 
J, K and B together determine the ground state prop- 
erties of the system instead of J and B in the isotropic 
Heisenberg model. This is a quite different feature from 
the previous studies. For example, we may choose K, 
J and B such that \^i) is the ground state of the sys- 
tem, it is entangled state but not a maximally entangled 
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FIG. 1: (Color on line) The thermal entanglement vs. the 
external magnetic field B and temperature. We choose the 
negativity as the entanglement measure, B and T are plotted 
in units of the Boltzmann's constant ks- At T = 0, N has a 
sharp transition from 1 to as B crosses the critical value of 
B = 3/2(J - K). The B = 0.35 line pointed out in the figure 
shows that for certain value of B it is possible to increase N by 
increasing temperature T. Figure l-(a), (b) are for different 
K, (a) K = -0.6 and (b) K = -0.7. 



one. The critical magnetic field can be increased by in- 
creasing the nonlinear coupling \K\, as shown in figure 
l-(b), where the entanglement is plotted as a function of 
T and B with the same parameters as in figure l-(a), but 
K = —0.7. It is obvious that the threshold temperature 
above which the entanglement vanishes has also been in- 
creased. More clearly, this point was shown in figure 2, 
where we plot the thermal entanglement in the system as 
a function of K and T. For B = 0, the thermal entan- 
glement is not zero only for \K\ > 0.4, as the figure 2-(a) 
shows, this indicates that \K\ > \ J\ is a necessary condi- 
tion for thermal entanglement to exist. Further analysis 
shows that for B = and with a specific temperature T 
the condition for thermal entanglement to exist is 



R < j- -l n ( _ ). 



(8) 



The condition Eq.© also holds for J > and K > 0, 
since we made no constraint on the derivation of Eq. (JSJ . 
On the other hand, it gives rise to a analytical expression 
for the threshold temperature T c in the case of B = by 



FIG. 2: (Color on line) This figure depicts the thermal entan- 
glement as a function of the temperature T and the nonlinear 
coupling constant K with a specific J = —0.4. The degree of 
entanglement is measured in terms of the negativity. Figure 
2-(a), (b) are for different B, (a) 5 = 0, and (b) B = 0.5. All 
parameters are re-scaled in units of fcg. 



solving K = J — -yln( 5+3e 2 1 " ), it shows again that T c 
depends on J and K. In addition to the above feature, 
there is a evidence that the threshold temperature is a 
monotonous function of \K\ with B — 0, the larger the 
nonlinear coupling constant \K\, the larger the thresh- 
old temperature. This point would be changed when the 
magnetic field is present (Fig. 2-(b)), as you see the 
threshold temperature is no longer a monotonous func- 
tion of \K\. We now look at the dependence of the ther- 
mal entanglement on J and K, with a fixed B and T, 
the selected results for this dependence were illustrated 
in figure 3. As figure 3 shows, the thermal entanglement 
may exist only in the regime of larger \K\, or smaller \J\. 

Now we discuss the experimental feasibility for observ- 
ing the thermal entanglement in the optical lattice. We 
may choose 23 Na trapped in an optical lattice as the 
system. Bose-Einstein condensates of unpolarized 23 Na 
have already been achieved by the MIT group [icj. The 
23 N a atoms have hyperfine spin 1, and the interaction 
among them is antiferromagnetic (t > 0) that is essential 
for the thermal entanglement to exist in isotropic Heisen- 
berg model. There are two parameters, the hopping ma- 
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FIG. 3: (Color on line) Thermal entanglement vs. the cou- 
pling constants J and K. The selected results are chosen 
from the regime of J, K < 0, which is relevant to atoms 23 N a 
in an optical lattice. The parameters chosen are B — 0.8 
and T = 0.2. All parameters are rescaled in units of fes, the 
Boltzmann's constant. 

trix t and the Hubbard repulsion U% (or Uq), we may 
control separately by adjusting the strength of the sinu- 
soidal potentials via the intensities of the laser beams, 
the coupling constants for the linear interaction J and 
the nonlinear term are hence adjustable. The two wells 
in the optical lattice in question can be isolated by in- 
creasing the barriers connecting to the other wells in the 
optical lattice. With these done, the two 23 Na atoms 
in the Mott regime then end up in a thermal entangled 
state. 

To conclude, we have studied the thermal entangle- 
ment in a optical lattice, the dependence of the ther- 
mal entanglement measured by the negativity on the lin- 
ear coupling constant J, the nonlinear coupling constant 
K, the external magnetic field and temperature was pre- 
sented and discussed. There are two different points in 
contract to the isotropic Heisenberg model, one is the di- 
mensionality and another is the nonlinearity of the cou- 
plings. When the nonlinear coupling is zero, there is no 
thermal entanglement in the regime J > 0, this is sim- 
ilar to that in the isotropic Heisenberg model in spite 
of different dimensionality. The nonlinear coupling re- 
ally change the feature of the thermal entanglement, in- 



creasing | if | increases the critical magnetic field and the 
threshold temperature. The thermal entanglement in an 
optical lattice with more than two coupled wells remains 
untouched. In future, we will investigate these problems 
and study how to map this natural entanglement onto 
photons and use it as a resource in QIP. 
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